Getting the elastic scattering length by observing inelastic collisions in ultracold 

metastable helium atoms 
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We report an experiment measuring simultaneously the temperature and the flux of ions produced 
by a cloud of triplet metastable helium atoms at the Bose-Einstein critical temperature. The onset 
of condensation is revealed by a sharp increase of the ion flux during evaporative cooling. Combining 
our measurements with previous measurements of ionization in a pure BEC, we extract an improved 
value of the scattering length a — 11.3 jl fg nm. The analysis includes corrections taking into account 
the effect of atomic interactions on the critical temperature, and thus an independent measurement 
of the scattering length would allow a new test of these calculations. 
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Understanding and testing the role of interparticle in- 
teractions in dilute Bose-Einstein condensates (BEC) is 
an exciting area of current research. Although measure- 
ments of the interaction energy and the spectrum of ex- 
citations of a BEC have confirmed the validity of the 
Gross-Pitaevskii equation |j| , there are still relatively few 
quantitative tests of other aspects, such as the effect of in- 
teractions on the value of the critical temperature (T c ) or 
the condensed fraction 0, El • The success in condensing 
metastable helium atoms (He*) 0,0, was greeted with 
interest in the community partly because the metastabil- 
ity offers new detection strategies unavailable with other 
species. To fully use these strategies however, we are still 
missing an accurate value of the s-wave scattering length 
a, the atomic parameter which determines all elastic scat- 
tering behavior at low energies. An accurate value of a 
would also be useful to help clarify some puzzling re- 
sults concerning measurements of He* in the hydrody- 
namic regime, in which two different ways of measuring 
the elastic scattering rate appeared to be in contradiction 
0. Also, because He is a relatively simple atom, theoret- 
icalpredictions of a are already in a rather narrow range 
0, and these calculations should be tested. 

A straightforward method to determine a is to use 
ballistic expansion of a BEC to measure the chemical 
potential for a known atom number. This was done in 
Refs. 0, 0|, but the measurements were limited by the 
calibration of the number of atoms, or equivalently the 
density of the sample. The reported values for a are 
20 ± 10 nm and 16 ± 8 nm, respectively. A more recent 
estimate, limited by similar effects, is a — 10 ± 5 nm @- 
In this paper we report a new measurement of a which 
makes extensive use of a unique feature of He*, sponta- 
neous Penning ionization within the sample. 

We exploit two specific situations in which the abso- 
lute atom number N is simply related to a and measured 
quantities: (i) for a pure BEC, the number is deduced di- 
rectly from the chemical potential fi and a, (ii) for a cloud 
at the Bose-Einstein threshold it is simply related to the 



critical temperature T c . Both /x and T c are accurately 
deduced from time of flight (TOF) measurements. Com- 
parison of ion rates from a pure BEC of known chemical 
potential and from a cloud at T c allows us to extract a 
and the ionization rate constants. The deduced value of 
a is independent of the absolute ion detection efficiency, 
assuming that this (poorly known) efficiency is the same 
in the two measurements. The ion signal is also used 
in another novel way: Since it provides a real-time ob- 
servation of the onset of BEC [l(| , we use it to reliably 
produce a cloud at the condensation threshold. 

A dense cloud of He* produces a steady flux of ions 
due to various ionization processes. Density losses due 
uniquely to ionizing collisions depend on the local density 
n according to: ( . . . — — — - 8 n 2 - L n 3 with n 

° V at I ionizing Ti ' 1 

the lifetime due to ionizing collisions with the background 
gas and (3 and L the 2-body and 3-body ionization rate 
constants defined for a thermal cloud. The total ion rate 
from a thermal cloud is given by: 
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The numerical factors reflect the fact that although 2 (3) 
atoms are lost in 2-body (3-body) collisions, only 1 ion is 
produced. Ionization measurements on a pure BEC were 
reported in |ll| . and, as a was not precisely known, 3 
and L were given in terms of a. 

For a precise measurement of a, corrections due to in- 
teractions must be taken into account. In the mean field 
approach, the density is given by [l|: 
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(2) 

with A(T) the thermal de Broglie wavelength, T the tem- 
perature of the cloud, fcs the Boltzmann constant, V the 
trapping potential energy, g = AirTi 2 a/m the interaction 

oc 

constant, fi the chemical potential and g a (x) = ar- 

i=i 1 

The ion rate at the phase transition <!> c can be derived 
from EclJEJ by a first order perturbation theory similar 
to Ref. [13] but with a fixed temperature rather than a 
fixed atom number. We use the chemical potential of a 
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gas in a harmonic potential at the BEC transition: 

/ic/fc S T c = |-+4 ff3/2 (l)f (3) 

This gives: 

$ c = (^) 3 x [+(1.20 + 2.48^+12.35^) 

+ ^(0.33 + 1.81^ + 6.75^) (4) 
+ ^(0.22 + 2.21^- + 6.50^) 

with ui = (2lo± +cl>||)/3, uj — (wiia;^) 1 / 3 , uj c — ksTc/h 
and A c = A(T C ). The numerical values come from the cal- 
culation of arithmetic series and are independent of any 
parameters of the cloud. The terms proportional to a/A c 
in Eq. @ account for the atomic interactions, while the 
corrections proportional to uj/uj c take into account the 
effect of finite sample size. For the typical parameters of 
our experiment we have (T c ~ 2 /^K and a = 12 nm) 
a/A c ~ uj/lu c ~ 0.02 corresponding to an interaction 
correction of 20%, 40% and 60% in the three successive 
terms in Eq. Q). Even though the first order corrections 
are large, we find, using an approach similar to Ref. |13| . 
that the second order corrections are negligible: —4%, 
1.8% and —3%, respectively, with the above parameters. 
In our case, finite size corrections are always smaller than 
those due to interactions; since a/A c ~ uj/lo c , the differ- 
ence corresponds simply to the numerical constants in 
Eq. £Q. 

Our setup has been described in Ref. jl 1| . Briefly, we 
trap up to 2 x 10 8 atoms at 1 mK in a Ioffe-Pritchard trap 
with a lifetime (r) of 70 s, and a lifetime due to ionizing 
collisions, (rj), estimated to be > 500 s. In a typical run, 
forced evaporation for 30 s cools a cloud to a tempera- 
ture near the phase transition. At this point, the rf-knife 
frequency is decreasing at a rate of 400 kHz/s. We use a 
clover leaf configuration with a bias field Bo = 300 mG. 
The axial and radial oscillation frequencies in the har- 
monic trapping potential are u!^/2w = 47 + 3 Hz and 
u)±/2tt — 1225 + 20 Hz respectively. A 2-stage, single an- 
ode microchannel plate detector (MCP) is placed 5 cm 
below the trapping region. Two grids above the MCP 
allow us either to repel positive ions and detect only the 
He* atoms, or to attract and detect positive io ns p ro- 
duced in the trapped cloud. As explained in Ref. [Tl| , to 
detect the ion flux, the MCP is used in counting mode, 
whereas we record the TOF signal at low gain (analog 
mode) to avoid saturation. As explained in Ref. [4(, the 
TOF signal is due to atoms in the m = state which 
are insensitive to the magnetic field. However, atoms in 
magnetic field sensitive states are still present, and their 
trajectories are affected by uncontrolled residual fields. 
Therefore, during the time of flight, we apply a mag- 
netic gradient in order to push these atoms away from 
the detector. The ratio between the detected atoms in 
the m = state and the initial number of atoms in the 
cloud is not well known so we use the TOF only to 
get the temperature. 
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FIG. 1: Variation of the ion rate as the atomic cloud is cooled 
through the phase transition for various initial densities (gray 
curves). The rf-knife frequency at t = is 2 MHz. The sud- 
den increase of the ion rate (crosses) occurs at the BEC tran- 
sition. The solid line passing through the transition points 
constitutes our empirical relation, named threshold curve. 



The crux of the experiment is to obtain a cloud of 
atoms at the phase transition. To identify the BEC 
threshold point, we monitor the ion signal. We have 
shown in Ref. [l(| that the onset of BEC is heralded by 
a sudden increase of the ionization rate associated with 
the increased density of the condensate. More precisely, 
the BEC threshold corresponds to the rapid change in 
slope of the ion rate vs time, or the maximum of the 2nd 
derivative [lj| . Figure Q shows a series of such ioniza- 
tion rates during evaporation through the BEC transi- 
tion. From these curves we can determine an empirical 
relation between the time of the onset of condensation 
and the ion rate preceding it. This relation stays valid 
only as long as we keep the same evaporation ramp and 
bias field. We will refer to this as the "threshold curve" 
in the following. Due to fluctuations of the bias field, we 
observe fluctuations of the time of BEC onset from run 
to run. These correspond to approximately ± 60 ms in 
time or ± 25 kHz in frequency, a value which agrees with 
independent measurements of the fluctuations of the bias 
field. 

Having established this relation, we can interrupt an 
evaporation sequence very close to the BEC threshold, 
and record the instantaneous ion rate as well as the cor- 
responding TOF signal. We typically chose the interrup- 
tion time for the evaporation ramp in advance, and then 
discarded the runs for which the interruption time was 
further than ± 60 ms away from the threshold curve. 

We fit the associated TOF spectrum to determine the 
temperature (Fig. We use Eq. J3J) together with /i c 
given in Eq. @ for the initial atomic density and then 
assume purely ballistic expansion of the cloud after the 
switching off of the trap. We shall refer to this fit as 
the Bose fit. The fits are weighted by an estimated un- 
certainty in each point of the TOF curve. To make this 
estimate, we chose a set of TOF spectra which appeared 
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FIG. 2: Time of flight signal corresponding to a cloud released 
from the trap (at t = 0) when its ion rate is on the threshold 
curve. We fit the data with an excluded window indicated by 
the vertical lines (width equal to the rms width of the TOF 
signal). A Gaussian function (dotted line) does not describe 
the central part of the data well, while the Bose function as 
defined in the text (solid line), does, indicating that the cloud 
is close to threshold. 



to show no systematic deviation from their fits and used 
them to estimate the amplitude of the noise. This noise 
varies as the square root of the amplitude of the signal 
indicating that we are limited by the shot noise of the 
atom detection. Our procedure is only an approximate 
indicator of the error bars. The chi square per degree of 
freedom \ 2 for the fits deduced in this way ranges from 
0.8 to 3. In addition, we first fit the TOF signal in its 
entirety and then exclude a successively larger window of 
the center of the spectrum out to the rms width of the 
spectrum. For all of our runs, we observe a variation of 
less than 5% and in most cases less than 3% of the fitted 
temperature as the excluded window is increased. 

In Fig. [21 we show an example of a typical TOF spec- 
trum and two fits to the data, using a Gaussian and the 
Bose function described above. Both fits shown use only 
the wings of the distribution. The ability of the Bose 
function to reproduce the center of the distribution with- 
out including it in the fit, unlike the Gaussian, confirms 
that the cloud is indeed close to the condensation thresh- 
old. In the following, we use the temperature given by 
the fit with an excluded window of half the rms width 
of the TOF signal in order to avoid the possibility of a 
small condensate component or other high density effects 
distorting our analysis. 

Before plotting the ion rate as a function of the critical 
temperature, we must correct the observed temperature 
to account for the hydrodynamic expansion of the cloud 
(see and references therein). This correction is done 
in the spirit of Ref. 15] which uses the Boltzmann equa- 
tion approach to take into account the effects of collisions 
during the expansion of the cloud. The collision rate in 
Ref. [13 is calculated using a Gaussian density profile. 
We rather use the value calculated for an ideal Bose gas 
|l(i| , which we have adapted to take interactions into ac- 
count. This correction depends on the scattering length 
but the effect on the final value of a is only of order 



0.3 nm for a ranging from 10 to 14 nm. We therefore 
simply assume a = 12 nm for this correction in the fol- 
lowing. In fact, due to the additional anisotropy of the 
expanding cloud in the horizontal (detector) plane, the 
fitting function should be modified; but a simulation of 
this effect shows that the correction to the temperature 
is less than 0.1 %. 

Finally, we correct the detected ion rate <E> Cj dct to 
account for the detection efficiency a such that <I> C = 
&c, det/o- It should be noted that the rate constants were 
obtained by ion rate measurements j 1 1 1 ] . This means that 
they were also corrected: j3 = fi&et/a' and L = Ldct/a'- 
Equation 2| shows that, as long as a = a', the detection 
efficiency cancels out and do not have any impact on the 
determination of a. We have checked experimentally that 
a = a' . To allow comparison with figures in earlier pub- 
lications, the ion rate and the theoretical curves shown 
in all the figures have been corrected using the same a 
as earlier, namely a = 0.42 [ToLITl). 

The results of this analysis are plotted in Fig. [3] 
Curves corresponding to the expected variation for three 
values of the scattering length are also shown. We see 
that a large fraction of the data falls between a = 10 
and 14 nm. The points at the highest temperatures how- 
ever, show a tendency to fall near the theoretical curve 
for a = 10 nm, while those at lower temperatures fall 
near a = 14 nm. To analyze this tendency further we 
examine the TOF fits more closely using the \ 2 value as 
an indicator of the confidence level of each measurement. 
A large \ 2 could mean that the Bose function with \i 
imposed to /i c is not the right fit function, and therefore 
that the cloud is not sufficiently close to T c . As shown 
in Fig. |3 outliers tend to be correlated with a large \ 2 ■ 
Note however, that the remaining scatter in the data is 
too large to be accounted for by our a priori estimates 
of the uncertainties in our ion rate or temperature mea- 
surements, and we presume that it is due to fluctuations 
in the determination of the BEC threshold. 

To determine the scattering length, we fit the black 
points in Fig. [3] with a as a free parameter and using (3 
and L parameterized by a as in Ref. 11]. The fit gives 
(all points are given equal weight) a = 11.3 nm. Our 
chief estimated uncertainty stems from the fact that our 
data show a systematic tendency to fall above the best 
fit at low temperature and below it at high temperature. 
To estimate this uncertainty, we fit the data (including 
gray points) separately for T c below and above 2 fiK. 
We find a = 13.8 nm for the low temperature data and 
a = 10.4 nm for the high temperature data. The uncer- 
tainties in the measurements of (3 and L also contribute 
to the uncertainty in Eq. 10} used for fitting. In fact, the 
uncertainties in (3 and L are highly correlated [Til ] and 
their contribution to the uncertainty is less than 0.5 nm. 

The error bars are obtained by summing quadratically 
the sources of uncertainties. Our final result for the scat- 
tering length is thus a = 11.3^'q nm. This result may 
be compared with the calculation of Ref. 8] . This work 
leads to a — 8 nm using the potential curves of Ref. ^} • 
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FIG. 3: Ion rate versus critical temperature. The points cor- 
respond to the results of 280 runs for which the ion rate was 
deemed sufficiently close to the condensation threshold. Gray 
indicates runs for which \ 2 in the TOF fits was above 2. The 
dashed line is the theoretical estimate for a — 10 nm, the dot- 
ted line for a — 14 nm (both including interaction corrections 
of Eq. [IJ. The two solid lines correspond to a = 12 nm, a) 
with interactions and b) without interactions, and illustrate 
the size of their effect. 

From Ref. [8j one also finds that a 0.5% shift of the re- 
pulsive part of that potential, would bring the theoretical 
value into agreement with our result. This 0.5% shift cor- 
respond to the estimated uncertainty in the potential of 
Ref. Another theoretical treatment gives a scat- 
tering length between 8 and 12 nm, also consistent with 
our results. 

Our result also allows one to give values for the 2 



and 3-body ionization rate constants. The error-bars 
of Ref. hj are modified to take into account the un- 
certainty of a. The uncertainty in the ion detection ef- 
ficiency also contributes to the uncertainty in the rate 
constants. As discussed in Ref. we will assume 

a = 0.42 to get the central values of the rate constants. 
We will include a one-sided contribution to the error- 
bars to account for the possibility, also discussed in [Tlj . 
that a could be a factor of 2 smaller. We finally get 
13 = 0.9±o;gXlO- 14 cm 3 /s and L — 2.5^ x 10~ 27 cm 6 /s. 
The rate constants are in good agreement with theoreti- 
cal predictions d, 0] . 

As shown in Fig. curves a) and b), our value of a 
is significantly shifted by the non-ideal gas corrections 
of Eq. @. Thus, when an independent measurement of 
the scattering length becomes available, our results can 
be used as a test of these corrections [R|. Note how- 
ever, that corrections to the critical temperature beyond 
mean-field theory [20j are small when one parametrizes 
the critical point in terms of an average density [lj. But 
an examination of the critical density measured in a lo- 
cal way, by imaging the ions from a cloud for example, is 
sensitive to critical fluctuation phenomena which go be- 
yond mean field theory similar to the homogenous case 
|20|. Thus, refinements of the ionization measurements 
described here promise to continue to provide new tests 
of BEC physics. 
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